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Abstract. The spontaneous symmetry breaking in a vibro-fluidized low-density 
granular gas in three connected compartments is investigated. When the total 
number of particles in the system becomes large enough, particles distribute them- 
selves unequally among the three compartments. Particles tend to concentrate in 
one of the compartments, the other two having the (relatively small) same av- 
erage number of particles. A hydrodynamical model that accurately predicts the 
bifurcation diagram of the system is presented. The theory can be easily extended 
to the case of an arbitrary number of connected compartments. 



1 Introduction 

A granular system is an assembly of macroscopic particles or grains interacting inelastically, 
i.e. mechanical energy is not conserved. These systems exhibit many peculiar behaviors as 
compared with molecular, elastic systems [T]. One of them is the tendency to spontaneously 
segregate into high and low density regions, i.e. to form density clusters |2|3|4|5| 6], Therefore, 
they provide prime examples of formation of structures in far from equilibrium systems. The 
physical origin of this effect is the inelastic character of the interactions between particles. The 
simplest form of clustering is presented by freely evolving granular systems, as a consequence 
of a long- wavelength hydrodynamic instability [2|3|4j . However this is an ideal theoretical state 
that can not be generated experimentally. 

Another interesting situation in which clustering effects show up, was observed several years 
ago in a seminal experiment 7]. A vertically vibrated system of grains was confined in a box 
separated into two equal connected parts by a wall of a certain height. For strong vibration, the 
particles distribute themselves on the average equally in the two compartments. Nevertheless, 
lowering the intensity of vibration (or the frequency) below a critical value, the spatial symmetry 
of the system is spontaneously broken, in the sense that a steady state is reached in which there is 
no equipartition of grains between the two compartments. Moreover, grains in the compartment 
with less number of particles have larger average kinetic energy than those in the high density 
part. By considering the exchange of particles between the two compartments as an effusion 
process, Eggers [8] proposed an analytical model to explain this phenomenon. Later on, the 
original experiment and also Egger's theory have been extended in different ways, considering 
systems with several compartments in which both clustering and declustering may occur [9], 
and also disperse mixtures of grains [10] . 

On the other hand, another symmetry breaking mechanism was reported in |llj . In this 
case, the equilibrium between the two compartments is hydrodynamic and not merely effusive. 
Consequently, this theory applies when the size of the opening connecting both sides is large as 
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compared with the mean free path of the granular gas in its neighborhood. This is the opposite 
limit of that for which a balance of the flux of particles is enough to guarantee stationarity. The 
extension to binary mixtures of this hydrodynamic equilibrium has been investigated in [12] , 
by means of molecular dynamics (MD). Let us also mention than while in Eggers' approach 
the external gravitations field plays an essential role, the mechanism put forward in ref. [llj is 
formulated in the zero gravitational field limit. 

The aim of this contribution is to extend the study carried out in [llj to the case of several, 
more than two, connected compartments. The size of the holes connecting the compartments 
is considered 'macroscopic', i. e. larger than the mean free path of the granular fluid next to it. 

It is a pleasure to dedicate this work to our good friend and colleague Carlos Perez. We 
know he would have enjoyed discussing this kind of structure formation and had apportted 
interesting suggestions as well as key ideas for its better understanding. 

2 The system and molecular dynamics simulation results 

The granular fluid is often modeled as a system of N smooth inelastic hard spheres (d = 3) 
or disks (d = 2) of mass m and diameter a. The inelasticity of the collisions is described by 
means of a constant, velocity independent, coefficient of normal restitution a, defined in the 
interval < a < 1. This rather simplified model will be the one considered here. It has proven 
to qualitatively capture many of the experimental features of real granular systems. Then, 
when two particles i and j collide having pre-collisional velocities Vi and Vj , respectively, their 
post-collisional velocities, v\ and v'j, are given by 

Vi=Vi — (er-wy)o-, (1) 

/ 1 + a 

v j = v o + ( - cr ■ v W tr > ( 2 ) 

where Wy = V{ — Vj is the relative velocity before collision, and cr is the unit vector directed 
along the the line joining the centers of the two particles at contact, away from particle j. 

The grains are enclosed in a box of width 35* and height L. For d = 2, S is a length, while 
for d = 3, it is an area. The box is divided into three equal compartments by two walls starting 
at a height h. A sketch of the system is given in Fig. [T] Collisions of particles with all the walls 
are elastic. To keep the system in a fluidized state, the wall located at the bottom is vibrated in 
a sawtooth way with a velocity Vb |13j . This means that all particles colliding with the wall find 
it moving 'upwards' with that velocity. Moreover, the amplitude of the vibration is considered 
much smaller than the mean free path of the gas in the vicinity of the wall, so the position of 
the wall can be taken as fixed. No external field acting on the particles is considered, something 
that can be understood as corresponding to the very strong shaking limit. 

We have performed a series of MD simulations of the system described above in two dimen- 
sions, i.e. for a system of hard disks. In all cases, the simulations started with the same number 
of particles in each of the compartments and a Gaussian velocity distribution. For given values 
of a and t>j, a steady state is reached in which, on the average, the particles are equally divided 
into the three compartments, as long as its total number iV is small. Nevertheless, when N is 
increased beyond a certain critical value, the spatial symmetry of the steady state is sponta- 
neously broken. More specifically, in all the simulations we have carried out, it is observed that, 
after a transient period, the particles concentrate in one of the compartments, the other two 
having roughly the same much smaller number of particles. Moreover, particles in the high den- 
sity compartment have significantly smaller velocity than particles in the low density ones. An 
example of a typical instantaneous snapshot is given in Fig. [21 which corresponds to a system 
of 500 inelastic hard disks with a = 0.9. The height of the system is L = 140er, and the width 
of each of the compartments S — 100cr/3. It is clearly seen that the number of particles in the 
central compartment is much larger than in the other two. The average density of the system 
is quite low. In particular, a dilute gas theory can be expected to apply in the lower part of the 
system, where there is no wall separating both compartments. 
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Fig. 1. Schematic picture of the set up considered. The wall at the bottom is vibrated in a sawtooth 
way with velocity i>b and small amplitude, while all the other walls are at rest. Collisions of the particles 
with all the walls are elastic 




Fig. 2. Instantaneous snapshot of a typical configuration of a system in the steady state with broken 
spatial symmetry. The values of the parameters are N = 500, a = 0.9, L = 140cr, 5* = 100a/3, h = 75a, 

and Vb = 2-\ 



' 2T(0) 



To have a clearer idea of what actually happens in the system, the time evolution along a 
trajectory of the population of each of the compartments for the same system as considered 
in Fig. [5] is plotted in Fig. [3J Time has been scaled with the inverse of the initial Boltzmann 
collision frequency vq(0) given by 

2Na /T(0)\ 1/2 

where T(0) is the (arbitrary) initial granular temperature of the particles. As usual, the granular 
temperature is defined from the average kinetic energy with the Boltzmann constant set equal 
to unity. It is observed that the symmetry of the system is broken very fast, and the system 
quickly evolves to a steady state in which the central compartment has about 400 particles, 
while the other other two compartments have about 50 each. It is worth to mention that in the 
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Fig. 3. Time evolution of the number of particles in each of the three compartments for the same 
system as considered in Fig. [2] Time is scaled with the inverse of the frequency given in Eq. (|3} 



figures shown here in relation with the time evolution of the system, only a relevant part of 
the latter is shown. The times reached in the simulations have always been much larger than 
those reported here in the figures, namely up to be really sure that the system was actually in 
a steady situation. 

As already mentioned, in all the carried out simulations, it has been found that when the 
spatial symmetry is broken, the particles agglomerate in one of the compartments, the other 
two having approximately the same (low) number of particles. Nevertheless, in some cases in 
which the symmetry breaking proceeds rather slowly, the final state is reached after going 
through an intermediate time-dependent unstable state. For a given coefficient of restitution, 
this usually happens when the number of particles is high enough. In the intermediate state, 
two compartments have the same relatively large number of particles, while the number of 
particles in the third one is significantly lower. After some time, this configuration decays, and 
the number of particles in one of the populated compartments decreases until reaching the level 
of the other low density compartment. An example of this kind of behavior is provided in Fig. 
|U The population of the central compartment grows for a while above the average, having the 
same value as that of the left compartment. Afterwards, the number of particles in it decreases 
until reaching the same value as in the right compartment. 

The symmetry or asymmetry of the steady state can be characterized by the set of param- 
eters 



with Ni being the average number of particles in the compartment i in the steady state. Of 
course, in a symmetric configuration, it is e; = 1/3 for all i. The resulting bifurcation diagram 
is shown in Fig. [5l where the three parameters are plotted as a function of a dimensionlcss 
control parameter £ m , that is proportional to the total number of particles and will be defined 
below. The data for different values of a are seen to collapse on the same curves. Moreover, it 
has been found that the bifurcation diagram is not altered by modifying the velocity of the wall 
Vb, as long as it is large enough as to keep all the granular system fluidized. For small values of 
£ m , the three compartments are equally populated, but for £ m > £o, one of the compartments 
has a larger number of particles than the other two. Note that the same symbol is used in the 
figure for the three compartments since they are equivalent, in the sense that they interchange 
their populations in different trajectories of the same system. 
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Fig. 4. Time evolution of the number of particles in each of the compartments for a system of TV = 900 
particles. The values of all the other parameters are the same as in Fig. [2] Note the slow relaxation of 
the system towards the steady state in this case as compared with that shown in Fig. [3] 




Fig. 5. Bifurcation diagram of the vibrated system with three compartments sketched in Fig. \T\ 
showing the relative average number of particles in each compartment = Ni/N, as a function of 
the dimensionless parameter £ m defined in the text. The symbols are MD results for the values of a 
indicated. The other parameters of the system are: L = 140<r, S = 100cr/3, and vt = 2-y/2T(0)/m. 
The lines are the theoretical prediction from the model developed in the text, where its meaning is 
explained 
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3 The model 

Consider the steady state reached by a vibrated dilute granular gas with only one compartment 
of width 3S and assume there are gradients only in the vertical direction, taken as the x axis. 
It is convenient to define a dimensionless space scale by 

where \{x) is the local mean free path given by 

\{x) = [C d a d - 1 n{x)Y l , (6) 
with C2 = 2v2, C3 = 7rv2) n(x) the local number density, and 

a[a> (d + 2fC 2 d n (d/2) [k* (a) - n* (a)} ' 1 ' 

Here n*{a) and n*{cx) are the dimensionless transport coefficients characterizing the Navier- 
Stokes heat flux in a granular gas |14j . Finally, C*( a ) is t ne dimensionless cooling rate due to 
the energy dissipation in collisions. The explicit expressions of these quantities can be found in 
ref. [T5]- For elastic systems with a — ► 1, k* tends to unity, while fi* and £* vanish. 

In the upper limit of the system x — L, £ takes its maximum value £ m = £(x = L) = 
■J a(a)Cd<J d ~ l N x , with N x = N/3S being the number of particles per unit of section of the 
vibrating wall. This is the quantity used in the horizontal axis of Fig. [5] In the Navier-Stokes 
approximation, the hydrodynamic pressure of the system is uniform and can be expressed as 

US 

V= — „ /(e m ), (8) 

where T is the temperature of the gas next to the vibrating wall and the function 

cosh ^ 

has been introduced. This function is plotted in Fig.[Sl It presents a maximum at £ = £oj where 
£0 ~ 1.20 is the non-zero solution of the equation £0 tanh£o = 1. For large values of £ it tends 
asymptotically to 2. 

Now come back to the steady state reached by the system with the three compartments 
represented in Fig.[TJ By extending the ideas developed in [11], we model this state by treating 
the compartments as independent, sharing only a thin, but macroscopic, layer of granular 
gas next to the vibrating wall at the temperature To. Macroscopic, hydrodynamic equilibrium 
between the three compartments requires that the pressure be the same in all them, and use of 
Eq. © gives 

/(d 13 ) = /(d 2) ) = /(d 33 ), (10) 

(i) 

where Qn is the maximum value of £ in the compartment i 1 i.e. next to the upper wall, 

^EE^RO^- 1 ^. (11) 

Since the total number of particles in the three compartments is N, the quantities must 
verify the condition 

4 X) + 4 2) + 4 3) = 3£ m = Sv^RC*^- 1 ^. (12) 

Of course, Eqs. (JTTJJ) always have the symmetric solution = = £m 3 = £m, in which the 
three compartments have the same average number of particles. From the observation of Fig. 
[H it follows that this is the only solution for £ m < £0 — 1-20. On the other hand, for £ m > £0 
other solutions are possible since there are two different values £~ and £ + , £~ < £ + , for which 
/(£;~) = /(£ + ). Therefore, in this parameter region two asymmetric solutions are possible: 
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Fig. 6. The function /(£) defined in Eq. It characterizes the dependence of the pressure on the 
total number of particles in the one-compartment system 

— Asymmetric steady state I. Two compartments have the same average number of particles, 
that is smaller than the average number of particles in the third one. Mathematically, it is 
defined by 

2C (/) +e+ (/) =U- (13) 

— Asymmetric steady state II. Two compartments have the same average number of particles, 
that is larger than the average number of particles in the third one. It is defined by 

C (//) +2a (//) =£ m . (14) 

Given the value of £ m > £o characterizing the system under consideration, the two above 
asymmetric solutions are found by numerically solving the equation /(£ _ ) = /(£ + ) together 
with Eq. (fl3|) or (fT4|) . respectively. This solutions are the lines plotted in Fig. The upper 
(lower) solid line is the function (£m^ 77 ^)- The dashed line is the function £m^ 77 ' > . The 

function ^m 11 ^ has not been plotted for the sake of clearness, taken into account that in the 
simulations, the system never was found in the asymmetric steady state II. Therefore, it seems 
that this state is unstable. It is seen that the agreement between the predicted values of the 
relative populations in the state II and the simulation results is pretty good. In particular, it 
is observed that their dependence on the coefficient of restitution a is scaled out when £ m is 
used as the control parameter. 

The asymmetry of the stable state / increases very fast as the total number of particles in the 
system increases. Even more, it can be observed in the figure that the steady average number of 
particles in the two less populated compartments decreases as more particles are added to the 
system. Let us stress that, even when the number of particles in these compartments is rather 
small, the hydrodynamic description presented here describes quite accurately the asymmetry 
of the system. 



4 Conclusions 

It has been shown that the problem of the spontaneous symmetry breaking in a vibrated 
granular gas with three compartments connected by a macroscopic hole, can be treated in a 
way similar to the case of two compartments. In fact, the theory can be easily extended to 
deal with an arbitrary number of compartments. The only important point to be verified is 
that the hydrodynamic pressure and temperature near the vibrating wall are the same in all 
the compartments. In the simulations presented here, it was found that they agree within the 
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numerical precision of the measurements. Of course, the exhibited symmetry breaking becomes 
much richer as the number of compartments is increased. 

An important point not yet well understood is why the system actually chooses the asym- 
metric state I with preference to the symmetric state and to the asymmetric state II, i.e. why 
the latter two are unstable. In principle, relevant information might be obtained by carrying out 
an stability analysis of the solutions. Nevertheless, such analysis is rather involved and has not 
been completed up to now. On the other hand, it is possible to compute the power dissipated 
in collisions for each of the states. This is done by taking into account that the energy balance 
requires that, in the steady state, the dissipated power be the same as the one injected through 
the vibrating wall. The latter for a system with a compartment is given by VbSp. For a system 
with three compartments, the powers dissipated in each of them have to be added up. Then, a 
simple analysis shows that it is maximum for the asymmetric state /. Also in the system with 
two compartments the asymmetric state dissipates more energy than the symmetric one. 

It is worth to mention that some care is needed when studying the state discussed in this 
paper both by numerical simulations and also by experiments. The steady state of the system 
with one compartment discussed at the beginning of Sec. [3] is known to be unstable when its 
width is larger than a certain critical value. More precisely, the system experiments a continuos 
spontaneous symmetry breaking in the direction perpendicular to the heat flux, next to the 
wall opposite to the vibrating one [16117] . The existence of this instability must be taken into 
account when deciding the sizes of the system to be used. 

This research was supported by the Ministerio de Education y Cienci'a (Spain) through Grant No. 
FIS2005-01398 (partially financed by FEDER funds). 
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